Application to the Binomial Summation of a Laplacian J 
Method for the Evaluation of Definite Integrals 

by e. c. molina 

Introduction 

THE numerical evaluation of the incomplete Binomial Summation, 
a problem of major importance for many statistical and engi- 
neering applications of the Theory of Probability, is a question for 
which a satisfactory solution has not as yet been obtained. Several 
approximation formulas have been presented, 2 each of which gives 
good results for some limited range of values of the variables involved; 
but a formula of wide applicability is still a desideratum. 

The purpose of this paper is to submit for consideration an approxi- 
mation formula which seems to meet the situation to a measurable 
extent. The writer derived it by applying to the equation 

f.T'-(l-.v)"-^ 

(i) E(*W -/>)--* =4 . 

Jo 

a method which is peculiarly efficacious for approximately evaluating 
definite integrals when the integrands contain factors raised to high 
powers. 

The method used constitutes the subject matter of Chapter I, 
Part II, Book I of Laplace's "Theorie Analytique des Probabilites." 
Poisson applied the method to the integrals in the equation 

z=n r *— /(i + x)^dx 

(2) £ r)p*(i - *>~ -™ 

x=eX ' x»-ia+xY+ l dx 



/•CO 

I * 

i/Q 



and published a first approximation, together with its derivation, in 
his "Recherches sur la Probabilite des Jugements." Poisson's ap- 
proximation seems never to have been used and was less fortunate 
than his famous limit to the binomial expansion which also was lost 
sight of until it reappeared under the caption "law of small numbers." 

1 Presented before International Congress of Mathematicians at Bologna, Italy 
in September, 1928. 

2 For an excellent resume of some well-known formulas, together with a discussion 
of their limitations, reference may be had to C. Jordan, "Statistique Mathematique," 
articles 37 and 38. 
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While the integrals in equations (1) and (2) are well known equiva- 
lent forms for the complete and incomplete Beta functions, the equa- 
tions themselves are not so familiar although one or the other will 
be found in Laplace, Poisson, Boole (Differential Equations) and at 
least two other places. 

Approximate Formula 

The approximate formula derived from equation (1) and submitted 
herewith for consideration is 

(3) 2(:)«»-*)"-A£'-w-^ : . 

where Si is the ith approximation to the infinite series 

J2R s T- l £l + (5 - 1)?\- 2 + (s - l)(s - 3)7^* ■ • •] 
1 + Z-Ra.Cl-3-5 ••• (25 - l)]2-« 

r, - rV2, 

(5) r- = (n - 1) log^ + (c - 1) log^-i + (» " c) log^, 

and a = np; T to be taken negative when a < (c — l)w/(» — 1). 
The first, second and third approximations to the infinite series 5 are 

c j? „ _Rx + RoT _ J ?x + -R 2 r + ■%,(! + T 2 ) 

where 



i? x = 4[(« - r ) - (c - 1)]/3V2(«- - 1)(» - c){c - 1), 
* 2 - (l/6)[l/(n - c) + l/(c - 1) - 13/(* - 1)1 
i? 3 = - (4/15)i?![i? 2 + 6/(» - 1)]. 

It will be noted that R», \R\\ and \R 3 \ are symmetric functions of 
(n — c) and (c — 1). 

For the limiting case (Poisson's Exponential Binomial Limit) 
where n = oo , £ = but np = a, we have 

r 2 = 1 + (c - 1) log (c - l)/c + (a - c), 
2?! = 4/3V2(c - 1), 
2? 2 = l/6(c - 1), 
R 3 = - (4/15)i? 1 i? 2 . 



BINOMINAL SUMMATION 101 

Numerical Results 

Since it is easy to compute the binomial summation directly when 
either c or n — c is small, the practical value of an approximate 
formula depends on its efficiency for large values of these quantities. 

The analysis given below under the heading "Derivation of the 
Approximate Formula" indicates that the successive 2? s 's in the 
series for S decrease when yc — 1 and yn — c increase. Therefore, 
when these two quantities are large, a few terms of the approximate 
formula (3) may be expected to give satisfactory results. As a 
matter of fact, the formula gives good results when -yc— 1 and 
\n — c are not large. To confirm this statement the Tables 3 given 
at the end of this paper are submitted. In the 4th column of each 
table are given 10 6 times the true values of 



p -|!(*W -*>)"-*• 



In the columns headed Ai, A« and A 3 are given 10 6 times the differences 
between the true values and those obtained by applying formula (3) 
with the first, second and third approximations to 5 respectively. 
Table I in Czuber's " Wahrscheinlichkeitsrechnung " was used for 
evaluating the probability integral in equation (3). 

The range of values of P covered by the tables is such that at the 
lower end of each section P > .0005 while at the upper end P < .9995, 
except where this latter condition would call for a value of c < 2. 
Of course, a larger or smaller range might have been given. The 
decision as to this question was based on the fact that several writers 
on the theory of statistics, when dealing with the normal law of 
errors, speak of an error exceeding 3 or 4 times the standard deviation 
as being a very improbable event. In order to keep the number of 
pages required for the tables within reasonable bounds computations 
were made only for even values of c. 

The values of a = np used are such that each of the values p = 1/2, 
p = 1/10 and p = 1/20 occurs twice; likewise each of the values 
n = 100, n = 50 and n = 30 occurs twice. 

A greater degree of accuracy than that indicated by the tables can, 
of course, be obtained by working out and using Ri, R- a ■ • • ; for this 
purpose, recourse should be had to equation (12) below and the 
details immediately following it. The only practical limitation to 
the use of formula (3) would appear to be the number of places given 

3 1 am greatly indebted to Miss Nelliemae Z. Pearson of the Department of 
Development and Research both for supervising the work of my computers and 
contributing personally several sections of the tables. 
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by the existing tables for the probability integral. However, this 
difficulty is encountered only when P, or (1 — P), is small, in which 
case T is large and the integral 



L 



T 

e~ t2 dt 



may be readily evaluated by computing the first few terms of the 
series 

[er"/2 rV^][l - Tr 2 + (l«3)rr« - (1-3 -5) TV* • • •]. 

where, as above, T\ = T->l2. 

When P is very small, the difference c — a = c — np is relatively 
large compared to a, and for this latter case recourse may be had to 
the approximate formula published by the writer in the American 
Mathematical Monthly for June, 1913. 

Derivation of the Approximate Formula 
Following Laplace closely, let us set 

(6) y(x) = Ye-», 

where Y = y(x ) is the maximum value of y(x). Then 

( 7) £* •'£*■'(*)*• 

the upper limit T being given by the equation 

(8) y(p) = y(x )e-r\ 
Assuming dx/dt expanded in powers of / so that 

(9) dx/dt = £*W 

s=0 

and setting R 8 = D B+ i/Di, equation (7) reduces to 
C ydx = YD^Rs f t'e-"dt. 

Jo s=0 •'— oo 

Our fundamental equation (1) may now be written 

2 R g C t*er*dt 

»=0 •/— oo 



— w ZRsf 

8=0 J—OC 



fe-i'dt 
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Integrating by parts and separating the terms involving fe'^dt 
from the terms containing e~ 11 , we obtain equations (3) and (4). 
To determine R s = D a+ \jD\, note that equation (6) gives / = (log Y 

— log y) 1 ' 2 and set v(x) = (* — .r )/(log Y — log y) 1/2 so that x may 
be written in the form 

X = Xq + V(x)t. 

This form for x gives the expansion (Lagrange's Theorem for the 
simple case where f(x) = x; see "Modern Analysis" by Whittaker 
and Watson) 

_ y- t* I d 3 ~ l V a \ 

x " fcjll\<te- 1 ) X=I0 ' 

Comparing this expansion for x with the previous expansion (9) for 

dxjdt, we obtain 

Di = v(x ) 
and 

Dx Ks \s K*)' *e" /*=* o' 

Up to this point no particular form has been attributed to the 
function y(x). From now on we deal with the function which consti- 
tutes the integrand of the integrals in equation (1). 

The function y(x) = :v c_1 (l — x) n ~ c gives the expansion (log Y 

— log y) = (x — xo) 2 [A + A\(x — x ) + A*(x — x ) 2 • • •]» where x 
= (c — l)/(w — 1) is the value of x for which y(x) is a maximum and 

, = 1 \ d«*Qog Y- logy) 1 

A ° (s + 2)\[ dx>+* \ x=x , 

or 

We are now prepared to evaluate R„. Set 

g = A + Ai(x - x Q ) + A 2 (x — X ) 2 • • • 
g a = d'g/dx'. 



and 
Then 



v = g- 112 , 
dv 2 . 

dx~ = ~ rgl ' 

?~ = (3/2)g" 7 ' 2 [(5/2)^ - g 2 gj, 
%=- 1r % \*£ - 9&gig + 12gx 3 ]. 
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Therefore, since g s — s \A s when x = Xo, 

Ri = - A<r*' 2 A lt 

R 2 = (3/2)A<r 3 l(S/4)AS - AoAsl, 

R* = - 2.4 - 9/2 [^3^o 2 - 3A«AxA, + 2AJJ 

Substituting for A , A u A« and A 3 the expressions derived by giving 5 
the values 0, 1, 2 and 3 respectively in equation (11), we obtain for 
Ri, i?o and R 3 the functions of n and c given on page 2. 

For values of 5 greater than 3 the direct evaluation of (i'V^/dx" 
by successive differentiation becomes very tedious. It will be found 
much more practical to use the following procedure, 4 where D is a 
symbol of operation, A = A and b = A\. 

/,7»«-(«+l)/2\ 

.4 -.«*. = (iA!)(^L_) 

[A A -(H-D/2 "1 T ri A -(H-D/2 I 

or 

The following equations give the details requisite for the formation 
of R, to R s inclusive; A a can be computed from equation (11). 

Db = A«, D°~b = A z , D 3 B = A 4 , D A b = Au 

D*b = A 9 , D 6 b = A 7 ,D 7 b = A it 

Db- = 2A } A 2 , 

D 2 b 2 = 2A X A Z + AS, 

DW = 2A 1 A i -\- 2A»A 3 , 

ZW = 2AiA 6 + 2A»A 4 +~A a \ 

DW = lAiAe +[2A«A- a + 2A S A 4 , 

D 6 b 2 = 2AiA 7 + 2A,A 6 + 2A 3 A- a + AS, 

DV = 3AMs, 

DW = SAMs + ZAxA-r, 

D*b* = 3 AM 4 + 6AiAzA 3 + A 2 3 , 

♦See DeMorgan's "Differential and Integral Calculus," 1842, page 328, art. 214. 
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DW = 3AM b + 6A1A2A4 + 3AiA 3 2 + 3AM 3 , 

DW = 3AMo + 6AiA%A t + 6AiA 3 A A + 3AMi + 34 2 .4 3 2 , 

Z)/,i = 4AM2, 

D 2 b A = 4^i 3 v4 3 + 6i, 2 i 2 2 , 

D 3 6 4 = 4i4i 3 i4 4 + l2AMtA* + 4^^.. 3 , 

D 4 6 4 = AAM- a + UAMiA 4 + 6.4 Ms 2 + \1A X A?A 3 + ^ 2 4 , 

D6 5 = 5i4i 4 i4 2 , Z>6 6 = 6i4i»i4j, 

D 2 b 5 = SAMz + IOAM2 2 , D n -b 6 = 6AMa + ISAM**, 

D 3 b 5 = 5AM 4 + 20AM*A a + 10,4 iM 2 3 , D6 7 = 7AM*. 

To illustrate the use of the procedure given above, let us evaluate 
R4. We have 

= -(S/2)i4 - 7/2 U4)+(l/2)(5/2)(7/2)i4o- 9/2 (2i4ii4 8 +i42 2 ) 

-(l/6)(5/2)(7/2)(9/2)^ - ll/2 (3^ 1 M 2 ) 

+ (l/24)(5/2)(7/2)(9/2)(ll/2M - 13 /M 1 4 
or 

P 4 = (5/2)^ - 6 C- ^o 3 ^4 + (7/2Mo 2 (iM 3 + ^ 2 2 /2) 

- (l/2)(7/2)(9/2MoiliMi 

+ (1/24) (7/2) (9/2) (11/2)^!^. 

Tables Indicating Degree of Accuracy Obtainable by 
Use of Formula (3) for Evaluating 



p-2Q*r-»- 



Pi = 1st approximation, Ai = (P - Pi) 10°, 

P 2 = 2d approximation, A 2 = (P - P 2 )10 6 , 

P 3 = 3d approximation, A 3 = (P - P 3 )10 6 , 
a = »/>, 

P 2 = (7/ - 1) log^y^- + (c - 1) log^-i + (n - c) logj^ 
7 ~ 4= I *-"<*/. 
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TABLE I 



7(10") 



7 > (lQ r ') 



A, 



a = 1.5, n = eo, p — 



2 


+ .3074653 


668154 


442174 


15991 


9518 


-1348 


4 


- .7612106 


140849 


65643 


10816 


1989 


30 


6 


-1.5874105 


12386 


4456 


1641 


303 


8 


8 


-2.2985028 


577 


170 


103 


20 






a = 1.5, n = 30,p = .05 



2 


+ .2865166 


657333 


446458 


21266 


17382 


-2083 


4 


- .8219430 


122536 


60772 


3684 


5617 


- 97 


6 


-1.6966449 


8211 


3282 


- 113 


914 


40 


8 


-2.4640017 


246 


85 


- 24 


50 


7 



TABLE II 



ZOO") 



PC lO") 







a 


= 5, n = oo, 


£=0 






2 


1.5461442 


985613 


959576 


-1681 


2590 


-798 


4 


.6837566 


833222 


734978 


-2036 


1897 


-138 


6 


.0000000 


500000 


384044 


2986 


1036 


- 4 


8 


- .5960752 


199621 


133376 


4219 


616 


17 


10 


-1.1358169 


54105 


31832 


2128 


286 


13 


12 


-1.6349406 


10385 


5452 


604 


123 


2 


14 


-2.1027717 


1471 


692 


107 


11 


- 5 


16 


-2.5454242 


159 


68 


14 


1 


- 1 







a = 


5, n = 100, p = .05 






2 


1.5596227 


986295 


962920 


- 373 


3331 


-732* 


4 


.6839234 


833281 


742162 


639 


3248 


-141 


6 


- .0162780 


490817 


384001 


2889 


2108 


- 38 


8 


- .6310024 


186097 


127961 


1989 


1431 


- 8 


10 


-1.1912234 


46029 


28188 


566 


703 


1 


12 


-1.7124507 


7723 


4274 


75 


208 


2 


14 


-2.2138799 


914 


463 


2 


38 





16 


-2.6715388 


79 


37 


- 1 


5 










a 


= 5, n = 50, p = .1 






2 


1.5742756 


987005 


966214 


830 


3962 


-724 


4 


.6844600 


833471 


749706 


3289 


4437 


-212 


6 


- .0335708 


481067 


383877 


2599 


3022 


- 62 


8 


- .6689826 


172053 


122145 


- 334 


2076 


31 


10 


-1.2524740 


38258 


24538 


- 789 


954 


51 


12 


-1.7994619 


5467 


3220 


- 278 


244 


24 


14 


-2.3191412 


520 


285 


- 48 


36 


6 


16 


-2.8175965 


34 


17 


- 5 


3 






I-P-Psl > \P -Pi\. 
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/(!()») 



P(10 8 ) 



a = 10, » = oo,/) = 



IP -Pal > \P - Pi 



TABLE IV 



Ai 



2 


2.5879363 


999874 


999499 


- 47 


67 


-37 


4 


1.8406742 


995381 


989662 


- 533 


275 


-53 


') 


1.2386541 


960089 


932912 


-1532 


518 


-50 


8 


.7094191 


842135 


779778 


-1585 


547 


-27 


10 


.2274981 


626172 


542069 


79 


425 


- 8 


1 > 


- .2200272 


377838 


303223 


1786 


322 


1 


14 


- .6408864 


182375 


135535 


2077 


239 


4 


16 


-1.0401811 


70640 


48740 


1374 


146 


4 


18 


-1.4215063 


22199 


14277 


628 


68 





20 


-1.7875189 


5737 


3454 


216 


25 


1 


22 


-2.1402533 


1236 


699 


58 


7 





24 


-2.4813121 


225 


119 


12 


1 


- 1 







a = 


10, n = 100, 


p-.l 






2 


2.6528972 


999912 


999679 


- 3 


71 


-24* 


4 


1.8917619 


996268 


992164 


- 15 


432 


-39* 


6 


1.2715533 


963931 


942424 


278 


1070 


-45 


8 


.7213308 


846163 


793949 


997 


1349 


-36 


10 


.2161911 


620099 


548710 


1213 


1179 


-19 


17. 


- .2564838 


358406 


296967 


503 


982 


- 2 


14 


- .7042404 


159638 


123877 


- 222 


771 


10 


16 


-1.1320595 


54691 


39891 


- 376 


470 


13 


18 


-1.5434535 


14526 


10007 


- 222 


206 


9 


20 


-1.9410214 


3025 


1979 


- 80 


66 


4 


22 


-2.3267578 


500 


312 


- 20 


15 


1 


24 


-2.7022383 


66 


40 


- 4 


2 






/(10») 



f(106) 



a = 15, n = «>,/> = 



4 


2.6780004 


999924 


999788 


- 18 


10 


- 4 


6 


2.1229551 


998660 


997207 


- 141 


54 


- 9 


8 


1.6324888 


989520 


981998 


— 525 


146 


- 14 


in 


1.1843012 


953019 


930147 


-1066 


242 


- 15 


12 


.7670044 


860975 


815249 


-1198 


274 


- 10 


14 


.3737500 


701445 


636783 


- 515 


245 


- 4 


16 


.0000000 


500000 


431911 


582 


203 





18 


- .3574541 


306598 


251141 


1311 


168 


2 


20 


- .7009899 


160758 


124781 


1351 


132 


2 


11 


-1.0324325 


72134 


53106 


961 


89 


2 


24 


-1.3532229 


27826 


19464 


523 


49 


1 


26 


-1.6645241 


9287 


6184 


228 


22 





28 


-1.9672925 


2700 


1715 


82 


8 


- 1 


30 


-2.2623270 


689 


418 


25 


2 
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TABLE IV— Continued 
a = 15, n = 30,p = .5 



6 


2.6019552 


999883 


999837 


52 


24 


- 12 


8 


2.0184138 


997845 


997388 


559 


152 


- 77 


10 


1.4626537 


980704 


978613 


2676 


411 


-239 


12 


.9237039 


904277 


899756 


5946 


488 


-343 


14 


.3942720 


711436 


707667 


5025 


227 


-205 


16 


- .1313195 


426335 


427768 


-1916 


- 73 


72 


18 


- .6581761 


175978 


180797 


-6392 


-382 


301 


20 


-1.1915875 


45979 


49369 


-4405 


-497 


316 


22 


-1.7378702 


6991 


8062 


-1346 


-278 


149 


24 


-2.3057782 


555 


715 


- 191 


- 69 


32 


26 


-2.9097701 


19 


30 


- 11 


- 7 


3 



TABLE V 



/(!()») 



P(10») 



a = 25, n = *> , p = 



10 


2.6086661 


999888 


999778 


- 11 


3 


- 1 


12 


2.2291734 


999191 


998583 


- 51 


11 


- 2 


14 


1.8705496 


995920 


993531 


- 159 


30 


- 4 


16 


1.5289263 


984700 


977705 


- 364 


59 


- 6 


18 


1.2015564 


955365 


939522 


- 617 


90 


- 7 


20 


.8863972 


894998 


866422 


- 765 


109 


- 7 


22 


.5818753 


794717 


752697 


- 651 


110 


- 7 


24 


.2867455 


657452 


606120 


- 253 


101 


- 6 


26 


.0000000 


500000 


447076 


268 


92 


- 3 


28 


- .2791919 


346482 


299814 


678 


84 


1 


30 


- .5515253 


217703 


182105 


837 


75 


2 


32 


- .8175896 


123790 


100070 


761 


62 


3 


34 


-1.0778902 


63708 


49782 


561 


45 


3 


36 


-1.3328647 


29718 


22460 


350 


29 


2 


38 


-1.5828952 


12593 


9212 


189 


17 


2 


40 


-1.8283181 


4860 


3445 


90 


9 


2 


42 


-2.0694313 


1713 


1178 


38 


4 


1 


44 


-2.3065005 


553 


370 


15 


2 


1 



= 25, n = 50, p = .5 



14 


2.3698187 


999598 


999531 


80 


16 


- 9 


16 


1.9447371 


997023 


996699 


404 


57 


- 32 


18 


1.5274793 


984621 


983580 


1329 


126 


- 74 


20 


1.1159208 


942735 


940539 


2844 


170 


-113 


22 


.7083182 


841759 


838881 


3763 


140 


-109 


24 


.3031406 


665931 


664094 


2415 


60 


- 56 


26 


- .1010188 


443200 


443862 


- 873 


- 20 


19 


28 


- .5055162 


237333 


239944 


-3426 


-102 


87 


30 


- .9117246 


98634 


101319 


-3499 


-166 


117 


32 


-1.3211006 


30859 


32454 


-2056 


-157 


95 


34 


-1.7352770 


7063 


7673 


- 774 


- 91 


50 


36 


-2.1561545 


1147 


1301 


- 191 


- 33 


17 


38 


-2.5860897 


128 


153 


- 31 


- 8 


4 



